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ABSTRACT 

This paper considers an unbounded eccentric centrifugal impeller with an 
infinite number of log spiral blades undergoing synchronous whirling in an incom- 
pressible fluid. The forces acting on it due to coriolis forces, centripetal 
forces, changes in linear momentum, changes in pressure due to rotating and 
changes in pressure due to changes in linear momentum are evaluated. 


NOMENCLATURE 



D efinition 

Vector acceleration 
Impeller perturbation 

Thickness of two-dimensional impeller 
Vector pressure forces 

Vector surface forces on control volume 
Vector force of the shaft on the fluid 

Change in force in i direction due to perturbation in j 
direction 

Principal stiffness 

Dimensionless principal stiffness, (K/pU^^R^) 
Cross-coupled stiffness 
Dimensionless cross-coupled stiffness 
Number of blades 


This work was supported in part by NASA Lewis Research Center, Contract No. 
NAG 3-180 and in part by the University of Virginia Rotating Machinery and Controls 
Industrial Research Program. 


323 



0 


0 ' 

R 

R 1 

r 2 


*b 

u 


u. 


Rl 

B 


V 

W E 

W, 


TOT 


£ 

9 


e 2 

0 b 

e" 


p 

U) 


Center of flow field 
Geometric center of impeller 
Dimensionless radii ratio, (R 2 /R]) 

Inner radius of impeller 
Outer radius of impeller 

Radius of impeller with respect to center of flow field 
Average radial inlet velocity in fixed coordinate frame 
Radial velocity at point B 

Relative velocity of rotating coordinate system at point B 

Analytic velocity expression for force equation 

Velocity in rotating coordinate system at point B 

Vector velocity at any point in rotating coordinate system 

Blade angle (tang = IL /wR ) 

Rl 1 

Dimensionless eccentricity 
Angle measured from x axis 
Entrance angle 

Exit angle 

Blade sweep angle (9 = cot B£n (R ? /R x ) ) 

Angle of fluid particle along blade 
Fluid density 

Angular speed of the impeller 
1. INTRODUCTION 


The analysis of turbomachinery requires a combination of solid mechanics and 
fluid dynamics. Vibrations of a centrifugal pump or compressor are controlled by 
the bearings, by the shaft geometry, bv fluid forces on the impeller, by seals, 
and other factors. Often the nature and magnitude of the forces and stiffnesses 
generated by bearings and shaft geometry are fairly well understood and can, in 
general, be accurately modeled. Those forces acting on the impeller can rarely 
be adequately determined as yet. These forces arise from the interaction of the 
impeller with the driven fluid, and are often called "aerodynamic forces," "aero- 
dynamic stiffnesses" or "aerodynamic cross-coupling". 


Much of the large body of work on incompressible flows in impellers assumes a 
centered impeller with varying degrees of complexity. A few examples include one- 
dimensional velocity vectors [1], two-dimensional considerations [2,3], and full 
three-dimensional analysis [4,5,6]. These are not discussed further. 


Calculation of aerodynamic forces was first 
circumferential variations of static pressure and 


reported by Alford [7], considering 
efficiency in axial compressors. 
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Work by Black [8], Barrett and Gunter [ 9 ], Lund [lo] , and Gunter, et al. [ll] » have 
demonstrated the importance of aerodynamic forces to the safe operation of turbo- 
machinery, but have not, in general, addressed the calculation of these forces, or 
their associated stiffnesses. Further, in the work by Barrett [_ 9 J and Gunter [_llj , 
these forces are modeled by their cross— coupled stiffness only; consideration is not 
given to the principal stiffness of the fluid- impeller interaction. 

Several recent papers have sought to more accurately quantify these forces. 
Colding— Jorgensen [^12j used a simple source and vortex flow model for the impeller 
along with a series of vortex sheets representing the pump volute. He then super- 
imposed a uniform flow over this model, without clear physical explanation, to deter- 
mine the forces on the volute, as a function of the volute spiral angle. The flow 
field emerging from the impeller is assumed to correspond with the typical one-dimen- 
sional analysis, except that its location is perturbed. The force calculated by this 
method is the total force acting on the entire impeller/volute model for uniform flow 
- not the force on the impeller. 


The work by Shoji and Ohashi [l3] considers an ideal impeller whirling about 
its geometric center. Blades are modeled by vortex sheets using unsteady potential 
theory with shockless entry to predict the flow field within the impeller. Forces on 
the blades are summed from the pressures calculated by an unsteady Bernoulli equation. 
The calculation is two-dimensional without volute or vaned diffuser effects but does 
allow for a whirl frequency other than the rotational speed. 


The papers by Thompsen [U,15] outline the assumptions and general theory con- 
tained in a proprietary program. The technique uses unsteady flow theory and deter- 
mines the stability of the rotor stage based on nonsynchronous fluid excitations. 

The proprietary nature of the program limits the discussion of the details of the 
procedure. An actuator disc method for calculating forces is currently being use ^ 
Chamieh et al. [l6]; however, the work is incomplete and detailed results are lacking 


A paper by Jenny and Wyssmann [l7] suggests the aerodynamic excitation is 
approximately two orders of magnitude smaller than other destabilizing mechanisms. 
Their simplified analysis, however, considers only the variation in radial clearance, 
while neglecting the change in momentum of the fluid. In a radial flow machine, t e 
change in force due to clearance variations should be negligible compared to changes 
in force due to momentum variations since the former are perpendicular to the main 

flow path. 


A paper by Washel and von Nimitz [l8] gives an empirical formula for aerodynamic 
cross coupling forces in compressor impellers. It is used for comparison to the 
theory developed in this work. 


The purpose of this paper [l9] is to consider the nature and magnitude of aero- 
dynamic cross- coup ling forces and stiffnesses as generated by the change in fluid 
momentum in an eccentric impeller. The approach is to perform a control volume 
analysis on a perturbed, infinite bladed impeller undergoing a steady orbit at the 
rotational speed of the shaft. This assumes incompressible, ideal flow withou 
volute or diffuser effects. It further assumes that the impeller consists of an 
infinite number of infinitely thin blades of simple geometry; the use of which allows 
an analytic expression for the aerodynamic forces to be obtained. Boundary layer 
development separation, and secondary flows, however, are not considered. e 
forces calculated are generated by a perturbation of the impeller about the flow 
center and thus are related to principal and cross-coupled stiffness terms. 
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separate analysis using a perturbation velocity would have to be performed in order 
to determine principal damping terms for the impeller. 


2. VELOCITIES IN AN ECCENTRIC IMPELLER 


Figure 1 shows a perturbed centrifugal impeller where 0^ 
of the impeller and 0 is the center of the steady flow field, 
flow enters through an inlet pipe centered at 0 and the radial 
versely with radius away from this point. The eccentricity a 
dimensionless eccentricity 


is the geometric center 
It is assumed that the 
velocity varies in- 
is small such that the 


e = 


a 

R 1 


is much less than unity, 
neglected. The impeller 
ic spiral blades with an 
0^, an inner radius R \ , 


For this analysis terms of order e^ and higher will be 
consists of an infinite number of infinitely thin logarithm- 
gle 6* The blades are backward curved with a blade angle 
and outer radius R 2 . 


The analysis is to be performed in a coordinate system rotating at rotor speed 
w a ^ out the center of the impeller. Only synchronous whirling about the flow center 
0 will be considered. The fluid is assumed to follow the shape of the blades through 
the impeller passage and to exit parallel to the blades. The particle path in the ' 
absolute reference frame is shown in Figure 2. The impeller is purelv two-dimen- 
siunal and ideal, incompressible flow is assumed. 


The radial inflow to a centered impeller in the fixed (non- rotat i ng) coordinate 
system centered at 0 is 


u = U. 


lL 

R1 r 


(1) 


where the average inlet velocity is 
■ 2tt 


U. 


R1 2tt 


1 

>7 / u(Ri,f 

J 0 


l)d0i 


Here is the angle of a fluid particle measured along r = R, . The impeller 
rotates with angular velocity oj. The blade angle is given bv 

, ft U Ri 
tan 3 = 

W R1 

From Fig. 1, the length R £ is 
^ = Rj + a cosGj 

Thus, along an eccentric impeller centered at 0; the inlet velocity is 
U R] 

Rl [ _ 1 

(Rl + a cosB}) ^Rl 1 4- r cos0j (2) 
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From the binomial approximation 


1 

1 + e cos0 j 


1 - e cos0 j 


which leads to 


u = U B - U R1 (1 - ecoaSj) 


( 3 ) 


This is the assumed form of the flow entering the eccentric impeller at a point B 
around the impeller . 

From the geometry of the eccentric impeller, this velocity can be resolved into 
components normal to and tangential to the impeller (shown in Fig. 3) 

% * U B C ° S6 ' U R1 (1 " E “ S V 

IT = -n sin6 - “U (1 -e cosBj) e sin0 : — "*U„ c sinO] 

BT B Ki 


for 6 < < 1. The relative velocity of the rotating impeller is 
V fi = ojRg = rnR^l +ecos0i) 

along the inner radius of the impeller. Resolving this velocity into normal and 
tangential components gives 

V BN = V B Sin<5 ~ a)R l £ sin0 l = w a sin0 l 

V BT = V B COS<$ = wR l U +ECOS0!) 

In the rotating coordinate system, the velocity of point B is W g = U fi - or in 
components 


W DM = U_ (1 - e cos0i - Ecot B sinOp) O) 

BN K1 

W BT = -u) Rl (1 + c cos9 2 + e tan B sin0 A ) (5) 

Since the tangential velocity does not enter the control volume, it will be 
neglected in this analysis. 


From continuity,- the radial velocity must vary inversely with r in the 
interior of the impeller. Let 0" denote the angle of the fluid particle as it 
travels along the log spiral blade relative to the x axis in the rotating coordinate 
system. The relative velocity of the fluid is then 


w M - u 


Rl 


Rl 

r 


(1 -e cos0' - e cot 3sin0') 


( 6 ) 


Let 02 denote the angle of the fluid particle at the exit of the impeller. Figure 
2 shows that, in the rotating coordinate system, the fluid which enters the impeller 
at angle 0j also leaves the impeller at angle 0j - 0^. Thus 
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e 2 - 0j- e. 


The relative velocity at the exit of the impeller is, at 6 
^1 

W N2 U R1 “ E cos0 2 - ecot 8 sin0 2 ) 


= 0 2 . 


For flow along a blade, the relative magnitude is 


W. 


TOT 


[V + v) 


given by 


o£ biades is ■«—*« fonow 


tan 0 


- W T = 


w. 


N 


Then the magnitude becomes 

W TOT ^ + ^ tan 2w h 


(7) 


where W is obtained from Eq. (6). While an exact representation of the variation 
of 0 with radius would be logarithmic, the- approximation of a linear variation will 
e used to facilitate the development of an analytic expression. It is 


0'- 0 


r ' Rl 
R 2 - Ri 6 b 


TOT 


.vc ~ sr;r - • 

= W TOT + 6)i - cos (0 + 

The general velocity expression becomes 

*TOT (r ’ 0) = U R1 T-^ 1 + (tan0)2j J/ 2T - e c 


- e cot 8 sin 


cos (& - — ^ 

\ R2 - Rl ) 

( 9 ~ R z ~- Riyj)(s in ( 0 + S)i - cos (9 + 8)j^j 


( 8 ) 
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3. FORCE ON THE WHIRLING IMPELLER 


The linear momentum equation for an accelerating control volume [20] is 
I [a + 2w x ? + wx(oo x r) + wxr]p dv 

J cv 

= I V pdv + / V p(V , dX) 

J CV •'CS 


F - 


where 

? = net surface forces acting on the control volume 

s 

a * rectilinear acceleration of the moving reference frame with respect to 
the fixed frame 

lo = angular velocity of the moving reference frame 
V = velocity of a particle in the moving reference frame 
to = angular acceleration of the moving reference frame 
r = position of a particle in the moving reference frame 

For a reference frame rotating with constant angular velocity, this expression 
reduces to 


F = 
s 


[2to x V + to x (u> x r)] pdv + j^cs 


V p(V . dA) 


The surface forces consist of two parts 


F » F + F 

s SHAFT p 


(9) 


where 


]£ = force of shaft on control volume 

SHAFT 

? 

P 


= pressure force acting on control volume 


Bernoulli 1 s equation for a rotating reference frame is 


P - | \Z x r\ 2 + | | V| 2 = CONST 


The pressure force on the control volume is given by 


F - - 
P 


/ 

J cs 


dX 
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or, from Bernoulli's equation 


-f I 7 ' 2 ] 


dA 


The final cxprcoaion for the force exerted by the 3haft (inpeller) on the 
fluid in the control volume is 


^SHAFT = / CP 2 « * V + plux (w x r )~] dv 

J cv^“ 


*/ 

Jcs 


V p V . dX + 


L (t | “ x?| 


- f |V|^|dA 


■} 


The desired force here is the force exerted by the fluid on the shaft, 
negative of the above expression. 


It is the 




ON 

SHAFT 


■d 


L 


p2o)xVd v - / p(ux (w x r ) dv 

CV j ^•'cv 


Coriolis Force 


J 


Centripital Force 


where 


-/ 5 


p $ . dX 


Change in Linear 
Momentum 


•/ 

J \ 


2 I a x r 2 dA 




Change in Pressure 
Due to Rotation 


^c.s 


f |V| 2 dA 




Change in Pressure Due 
to Change in Linear 
Momentum 


r = (r + a cos0) (cos0 i + sin6 j) 

-*■ -*■ 
u) = wk 

V = W = U„, fl + (tanB) 2 /z 1 - e cos0' 

TOT R1 r j J L -i 

- ecot 6 sin0' Jjsin + ® )i - cos + I 
dA = (cos0 i + sin0 j) rd 0 dz 


($ . d 1) = W N - U ^ (1 - ecosB' - e cot Bsin0') 


(10) 
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dv » r d 0 dr dz 

Integration is taken over the control volume from to R 2 , and from 0 to 2n. On 
the control surface, integration is performed both at r = R 2 and at r = Rj over 0 
from 0 to 2 tt . The thickness of the impeller is assumed to be the constant b. 


For convenience in evaluation, the force acting on the shaft (impeller) 
due to the fluid is evaluated as five parts. 


£ = ~ CD - © ~ (3) -(4) +(5) 

ON cv ^cv cs ^-'cs 

SHAFT 


cs 


(id 


After integration, they are 
Coriolis Force 

= 2it aco 2 p b R 

* O' 


r*R - r "i 


1 2 1 

tan g 

I 0 , 

b -J 



sin (0, + 6) 4- sing - cotg cos (0, + B) + cotg cosBJi 

D D 


+ [cos (0, + 0) - cosg - cotg sin (0, + 0) 4- cotB sin0]j 
b d 


) 


Centripetal Force 
© = -y ira u' 2 pb (R 2 2 - Ri 2 )t 


Change in Linear Momentum 

2irao) 2 pb Rising (tang) 



[1 + (tan g) 2 ] 

cotg cos (0. + g) + sing 
R 2 b 

(0 b + g) - cotg sin (9 b + g) 
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Change In Pressure Due to Rotation 


(4) ■ 7T au 2 pb(R 2 2 - R x 2 )l 

Change in Pressure Due to Change in Linear Mom en t uni 

© =rr auj 2 pb R| 2 (tanB) 2 [1 + (tanB) 2 ] 



The physical significance is indicated above each term. 

The five parts can be summed together yielding the vector force acting 
on the shaft. Dividing the force into components gives F and F . A similar 
impeller eccentricity along the y axis would yield the forces F 7 and F . What 
is of interest here is the stiffness obtained from xy yy 


/¥ 



The principal stiffness is 


K xx = " ™ 2pb R l 2 ( 0 + 


+ 


m + Q- 


02 ) 


where 


m - 


-2 


1 


(R - 1) tanB [1 + (tan B) 2 r 2 


V? 


x [-sin (0^ + B) + sinB + cotB cos (6^ + 8) + cotB cosB] 

0 = \ (R 2 - i) 

QJ = - 2 sinB (tanB) 2 [1 + (tanB) 2 ] 

I 

x f~ = sin ( 0 5 + | cotB cos (0 fe + B) + sinB + cotB cosB 

E] = - (R 2 - 1) 

0 = (tang) 2 [l+ (tang) 2 ] 
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cos e, + — cote sin 0. + 1 

bp b 


The cross-coupled stiffness is 


K = - iroj 2 pb R± 2 ( A + + ^ ) 

yx 


where 


A = - 2 (R - 1) tang [1 + (tang) 2 ] ' 2 


[cos (0, +0) - cot0 sin (0, + 0)] 


/$\ = - 2 sin0 (tan 0) 2 [1 + (tan0) 2 ] 

x [- cos (0 + e) - * cote sin (6, + 6)] 


A = (tan B) 2 [1 + (tan B) 2 ] 

x [- “ sin 0^ — ~ cotB cos 0^ + cotB] 

Note that both principal and cross-coupled stiffnesses are obtained. In a rotor 
dynamics analysis, the other stiffnesses would be given by 

K = K 
yy xx 

K = -K 
xy yx 

For comparison with other works and convenience in plotting results, the 
stiffnesses have been made dimensionless in the following manner 

K ij = K ij 

K ij ” pU 2 R x puj 2 R 3 tan 2 B 

Ki 


The resulting expressions are 


— i isk +A +A +A +A) 


K yx " _7T R x tan 2 B~ + ^ 

The dimenionless stiffnesses are functions of the ratio b/Rj , the_radius ratio R, 
the blade angle B, and angular sweep of the blade 0 fe . The terms R, B, and 6 b are 
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related to one another so actually only two are independent. Also the terra b/R, 
appears only as a multiplicat ive constant. Thus only the radius ratio R and the 
blade angle 6 need be varied for these stiffnesses. 

4. RESULTS 


Figure 4 shows a plot of the dimensionless stiffnesses K and K vs. blade 
angle B for a radius ratio R = 1.5. For plotting purposes b/ft* has been set to 
unity. It shows that the blade angle has a strong influence on both the magnitude 
and sign of the stiffnesses. The principal stiffness is large and positive for 
low blade angles. Pumps generally have a low blade angle so they fit in this 
category. The cross-coupled stiffness is negative for this level of blade angle. 
The large positive principal stiffness is generally a stabilizing effect while 
C e " e S atlve K yy term is destabilizing. In the range of pump blade angles, 10 
to 20 degrees, the principal stiffness term is larger than the cross-coupled 
stiffness probably leading to an overall stabilizing effect. Compressors generally 
have large blade angles (as defined in this paper), perhaps in the range of 60 to 
0 degrees. With a radius ratio R = 1.5, both the principal stiffness and cross- 
coupled stiffness are negative indicating destabilizing effects. Field results 
indicate that pumps are usually stable while compressors are sometimes unstable. 
Thus the theory developed here appears to agree, at least qualitatively with 
results from actual machines. 


A se mi -empirical formulation for the cross— coupled stiffness was developed in 
[18] of the form 

6300 HP (Mol Wt)p Q 

K yx ~ NDhp s 

where 

HP =* pump horsepower 

Mol W£ ** fluid molecular weight 

P Q - discharge density 

p =* suction density 
s 

N = speed in RPM 

D = impeller diameter 

h = restrictive dimension in flowpath 

While this expression was developed for compressible flow, it can be reduced to incom 
pressible flow. Non-dimensionalizing this reduced expression gives 

K =* -0.6r (b/R ) n cot$ R(R-l) 

y* 1 

where 

n = the number of blades 

HP * (<i)R ) 2 - iuU_ R. cot 6 
2 R 1 1 
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The curve for this expression was added to Figure A with a blade number of eight 

in this paper agrees fairly well with the semi-empirical formula. 

Fieure 5 gives the results for a radius ratio R = 2.0. Again the principal 
stiffness is large for small blade angles. At larger blade angles, it decays 
t“"ro buJ doe/not go negative as it did for R - 1.5. The 

iar ger"nega t ive coefficient than the 

theory. 


One of the purposes of this study is to examine the effects of various terms^ 

, .iff ___ Table 1 gives the five terms involved in the pnnc p 

on the stiffnesses. Table 1 gives tn ^ angles, the dominant effect is 

ness vs. blade angle for R .5. Thig term arises in Bernoulli’s equation due 

the change in pressure ue larger blade angles, the dominant terms 

to the rotating coordinate system. A other term from Bernoulli’s equation is 

hprome the linear momentum ones. Now tne ocner uim ^ Tn} . 1p 2 

the largest but it is nearly balanced by the change in linear momentum. Table 2 
shows the three terms in the cross-coupled stiffness vs. blade angle for R 1. • 

A low blade angles, the Coriolis force produces the cross-coupled stif ness effect. 
At higher blade angles, the linear momentum terms become dominant, 
give the numerical values for R - 2.0. 


It is sometimes difficult to obtain a physical feel for the results of an 
analysis when dimensionless parameters are used. In an effort to impart a better 
idea of the numerical values developed from the theory here, some sample impellers 
were chosen. The results are given in Table 5. Generally small pumps will Mve 
very small stiffnesses while large pumps such as boiler feed pumps can have lair v 
large stiffness acting on them. The stiffness value here K xx - 124,000 lbf/m is 
the same order of magnitude expected from the seals in a pump. A small compressor, 
even at high speed, produces a relatively small cross-coupling stiffness. Large 
compressors have somewhat larger coefficients. Often such a compressor will have 
eight or ten stages with a fairly flexible shaft. The cumulative effect can produce 

an instability. 


5. CONCLUSIONS 


The theoretical solution for a simple impeller model shows the principal an 
cross-coupled stiffnesses to be of about the same order of magnitude. For blade 
angles (3) less than thirty degrees, common for water pumps, the flow see ®s t 
provide tabilization for the shaft, while for the larger blade angles, common 
for compressors, the impeller is generally destabilizing. Sine. I the ° £ 

if and K are nearly the same, it would be important to include both terms in 
any rotordynamics analysis of the shaft, rather than incorporating just the cross- 

coupling terms. 

• Though this work does not fully resolve the nature of aerodynamic forces 
centrifugal machines, it does offer more understanding of these forces ^ m P°* tant 
element of any future work, however, is the availability of experimental data for 
verification. As of the present, no such data exists, although efforts are being 
perused in this area. 
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TABLE 1 


EFFECT OF VARIOUS TERMS ON PRINCIPAL STIFFNESS FOR R = 1.5 


Blade 

Angle 

Coriolis 

Force 

Centripetal 

Force 

Change In 
Linear 
Momentum 

Change In 
Pressure Due 
To Rotation 

Change In 
Pressure Due 
To Change In 
Linear Momentum 

Principal 

Stiffness 

B 

0 

0 

0 

E 

0 

K 

XX 

10.00 

15.00 

20.00 

25.00 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

60.00 

65.00 

70.00 

75.00 

80.00 

-.75 
-.67 
-.57 
-.50 
-.45 
-.42 
-.40 
-.40 
-.41 
-.43 
-.47 
-.53 
- -63 
-.81 
-1.18 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

.63 

-.09 
-.15 
-.21 
-.30 
-.44 
-.65 
-.98 
-1.55 
-2.55 
-4.46 
-8.44 
-17- 7B 
-43.96 ! 

-140.25 ! 

-714.71 

A 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

-1.25 

.03 
.08 
.15 
.26 
.44 
.73 
1.20 
2.00 
3.44 
6.20 
12.00 
25.75 
64. 53 
207.92 
1066.62 

144.75 
59. 6‘/ 
29.76 
16.74 
10.05 
6.15 
3.62 
1.80 
.33 
-1.05 
-2.59 
-4.66 
-8.04 
-14. 94 
-34.20 


TABLE 2. - EFFECT OF VARIOUS TERMS ON CROSS-COUPLED STIFFNESS FOR R = 1.5 


Blade 

Angle 

Coriol is 
Force 

Change In 
Linear 
Momentum 

Change In 
Pressure Due 
To Change In 
Linear Momentum 

Cross- 
Coupled 
Stif fness 

Cross- 

Coupled 

Stiffness 

8 

A 

A 

A 








K 

xy 

Wachel [18] 

10.00 

15.00 

.33 

.71 

.03 

.10 


.06 

-42.96 

-64.14 

20.00 

.91 

. 18 


09 

-39.32 

-42.21 

25.00 

1.07 

.27 


12 

-28.73 

-31.07 

30.00 

1.23 

. 38 


15 

-21.52 

-24.25 

35.00 

1.41 

. 53 


18 

-16.84 

-19.59 

40.00 

1.64 

. 74 


19 

-13.69 

-16.15 

45.00 

1.95 

1. 05 


19 

-11.46 

-13.48 

50.00 

2.37 

1.53 


11 

-9.78 

-11.31 

55.00 

3.00 

2. 32 


11 

-8.40 

-9.49 | 

60.00 

65.00 

3.96 
5. 57 

3.71 
6. 45 

-2 

70 

31 

-7.10 

-5.62 

-7.92 ! 

-6. 53 

70.00 

75.00 

80.00 

8.52 

14.90 

33.13 

12.65 

30.06 

101.59 

-6 

-22 

-87 

-533 

93 

29 

79 

98 

-3.47 

.47 

9.66 

39.00 

-5.27 

-4.12 

-3.03 

-1.99 
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TABLE 3 


EFFECT OF VARIOUS TERMS ON PRINCIPAL STIFFNESS FOR R * 2.0 


Blade 

Angle 

Coriolis 

Force 

Centripetal 

Force 

Change In 
Linear 
Momentum 

Change In 
Pressure Due 
To Rotating 

Change In 
Pressure Due 
To Change In 
Linear Momentum 

Principal 

Stiffness 

8 

0 

0 

m 

E 

0 

K 

XX 

10.00 

-.89 

1. 50 

-.09 

-3.00 

.03 

247.42 

15.00 

-1.53 

1.50 

-.22 

-3.00 

.08 

138.95 

20.00 

-1.58 

1.50 

-.35 

-3.00 

.15 

77.74 

25.00 

-1.50 

1.50 

-.51 

-3.00 

.26 

46.85 

30.00 

-1.42 

1.50 

-.73 

-3.00 

.44 

30.17 

35.00 

-1.36 

1.50 

-1.06 

-3.00 

.73 

20.43 

40.00 

-1.33 

1.50 

-1.59 

-3.00 

1.20 

14. 35 

45.00 

-1.33 

1.50 

-2.46 

-3.00 

2.00 

10.35 

50.00 

-1.37 

1.50 

-4.00 

-3.00 

3.44 

7 . 60 

55.00 

-1.45 

1.50 

-6.92 

-3.00 

6.20 

5. 64 

60.00 

-1.58 

1.50 

-12.95 

-3.00 

12-00 

4.22 

65.00 

-1.79 

1. 50 

-27.08 

! -3.00 

25.75 

3.16 

70.00 

-2.15 

1.50 

-66.57 

| -3.00 

64.53 

2.37 

75.00 

-2.76 

1.50 

-211. 50 

I -3.00 

207.92 

1.77 

60.00 

-4.05 

1.50 

-1074,58 

| -3.00 

j 

1066162 

1.32 


TABLE 4. 


EFFECT OF 


VARIOUS 


TERMS ON CROSS-COUPLED 


STIFFNESS FOR R 


2.0 


Blade 

Angle 

Coriolis 

Force 

Change in 
Linear 
Momentum 

Change In 
Pressure Due 
To Change In 
Linear Momentum 

Cross- 

Coupled 

Stiffness 

Cross- 

Coupled 

Stiffness 

8 

A 

A 

A 

K 

xy 

Wachel [18] 

10.00 

15.00 

20.00 

25.00 

30.00 

35.00 

40.00 

45.00 

50.00 

55.00 

60.00 

65.00 

70.00 

75.00 

80.00 

-.37 

.44 

1.12 

I. 63 

2.07 1 

2.52 

3.03 

3.69 

4.57 

5.84 

7.79 

II. 00 
16.92 
29.68 
66. 16 

-.02 
.04 
.14 
.26 
.41 
. 61 
.88 
1.28 
1.89 
2-89 
4.68 
8.18 
16.11 
38.40 
130.04 

.[ 

.09 
.14 
.19 
.25 
.33 
.41 
.50 
. 59 
. 64 
.56 
.00 
-2.20 
-10.85 
-51.91 
-353.46 

_ 

30.89 

-26.93 

-34.58 

-31.07 

-26.49 

-22.65 

-19.69 

-17.44 

-15.71 

-14.31 

-13.05 

-11.60 

-9.23 

-3.65 

15.36 

-171.04 

-112.56 

-82.86 

-64.68 

-52.24 

-43.07 

-35.94 

-30.16 

-25.31 

-21.12 

-17.41 

-14.06 

-10.98 

-8.08 

1 
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Impeller 

Small Pump 
Large Pump 
Small Compressor 


Large Compressor 








F 1 gure 2 . 
frame [ 



- Fluid path through impeller over passage time (T). Fixed coordinate 
Ip = 10 T - 0 b]* 



Figure 3. - Velocity components at impeller inlet. 


3 



(b). Relative velocity between fluid and impeller. 



(c). Relative velocity components in rotating coordinate system. 


Figure 3. - Concluded. 
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